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Thermodynamics of ideal Fermi gas trapped in an external generic power law potential JJ = Tf i Cj| —| n< 
are investigated systematically from the grand thermodynamic potential in d dimensional space. These proper¬ 
ties are explored deeply in the degenerate limit (/i >> KbT), where the thermodynamic properties are greatly 
dominated by Pauli exclusion principle. Pressure and energy along with the isothermal compressibility is non 
zero at T = OK, denoting trapped Fermi system is quite live even at absolute zero temperature. The nonzero 
value of compressibilty denotes zero point pressure is not a constant but depends on volume. 

PACS numbers: 03.75.Ss, 05.30.-d, 05.30.Fk 


1. Introduction 

The constrained role of external potential can change the characteristic of quantum gases[H [2] [3] 0] [5] [6]. 
An increasing attraction towards this subject was noticed, after it was possible to create Bose-Einstein con¬ 
densation (BEC) in magnetically trapped alkali gases[7J[8j 9j. A lot of studies are performed to understand 
the behaviour of ideal Bose gas l(.) f 11, 321 [131 H] as we ll as ideal Fermi gasjlll, 12]- But unlike the Bose 
gas, Fermi gas does not condensate due to the Pauli exclusion principle, the question of a large number of 
particles occupying a single energy state does not even arise in this case. At sufficiently low temperatures, 
Fermi gas displays its own brand of interesting behavior{1T[ 32] . as its fugacity zf can take on unrestricted 
values: 0 < zf < oo m unlike the Bose gas who has restricted value of fugacity 0 < zb ^ m■ The 
behavior of thermodynamic quantities of Fermi gas are remarkably governed by Pauli exclusion principle. 
For instance, the ground state pressure also known as degeneracy pressurefll] hi ideal Fermi gas is nonzero 
unlike Bose gas or classical gaspTj. Nevertheless, a lot of efforts are made to understand different properties 
of Fermi system such as magnetism[13], conductivity[F6]. transport properties[T7], equivalence with ideal 
Bose gas Q21Q1, dimensionality effects [2D], degeneracy [21]. polylogarithmis [122]. q-deforemed svetm [251 [23] , 

In a real systems, of course interaction between particles does exist. But taking it into account makes 
the problem difficult to solve analytically. Neveretheless, to understand the effect of interactions in quantum 
gases and to retain the essential physics, we approximately represent the real system by non interacting par¬ 
ticles in the presence of an external potential [2l 1251 [26] . The trapping potential in atomic gases provide the 
opportunity to manipulate the quantum statistical effects. Some drastic changes are noted in case of Bose 
system]]], [2j under trapping potential. For instance BEC is possible in d = 2 in trapped Bose gas which was 
not a case in ideal Bose gasp], [2]. Therefore, it will be interesting to investigate how the trapping potential 
do change the properties of Fermi gas. Li et al. [?], in their work has presented internal energy, heat capacity, 
ground state energy of Fermi gas under spherically symmetric potential (U = br l ) in arbitrary dimension. 
However, in the present study we have investigated properties of ideal Fermi gas under generic power law 
potential U = X^=i hi d-dimension, which will be symmetric under certain choice n,, Xi and Cj. So, 

we can recounstruct the results of Li et al. [?] choosing this special condition. At first, we have calculated 
the density of states which enables us to determine the grand potential. Then from the grand potential of 
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the system, we have derived the thermodynamic quantities such as internal energy E entropy S, pressure 
P , number of particle N, helmholtz free energy A , isothermal compressibility kt, specific heat at constant 
volume Cy and pressure Cp and their ratio. In the high temperature limit, the thermodynamic quantities 
of free quantum gases reduce to the form of classical gas [11 J . Same trend is observed in case of trapped 
systempQ too. Therefore, the low temperature limit of the thermodynamic quantities of quantum gases are 
particularly important, as the true quantum nature is explicit in that region. In Bose system, low tempera¬ 
ture limit refers to condensed phase. It was found trapping potential changes the general criterion of BEC 
as well as the condition of jump of specific specific heatp], 2|. So, low temperature limit of thermodynamic 
quantities related to trapped Fermi system have been investigated using Sommerfield expansion) 11], It will 
be very intriguing to investigate energy and pressure of trapped Fermi gas while the gas is in the degener¬ 
ate limit to check whether they remain nonzero under generic trapping potential at T = OK. Isothermal 

compressibility (inverse of bulk modulus) is also calculated to check whether the ground state pressure has 

2 

volume dependency or it is marely just a constant. A point to note that, in the hamiltonian instead of 4i_ 
type kinetic part, we have took ap s , where p is momentum, a is constant and s is a arbitrary kinematic 
parameter. Different kinematic characteristics of quantum systems lead to different characteristics[T1 1281129j . 
So, choosing arbitrary kinematic parameter, we make the significant conclusions in more generalized way. 
In the current study it is found the concept of effective volume plays an important role in trapped Fermi 
gas, as seen in trapped Bose gas [I]. 

The report is organized in the following way. The density of states and grand potential are calculated 
in section 2. Section 3 is devoted to investigate the thermodynamic quatities. Properties of degenerate 
Fermi gas is presented in section 4. Results and discussions are presented in section 5. The report is 
concluded in section 6. 


2. Density of States and grand potential of Fermi gas under generic power law potential in d 

dimension 


Considering the ideal Fermi gas with kinematic parameter l in a confining external potential in a d- 
dimensional space with energy spectrum, 


d 

e(p,Xi) = bp l + ^2ci\ — \ ni (1) 

i= 1 (Ll 

Where, b, l, ai , c,;, n* are all postive constants, p is the momentum and x t is the i th component of coordinate 
of a particle. Here, c t . ai and rii determines the depth and confinement power of the potential. With l = 2, 
b = one can get the energy spectrum of the hamiltonian used in the literatures HD Eg EH E| . For the 
free system all n t —> oo. 


Density of states can be obtained from the following formula, 

d d rd d p 


P(t) = 




5(e - e(p,r )) 


So, from the above equation density of states is, 


where, 


p(e) = B 


r (f + 1 ) c x-i 

r(x) 


B = 


gVdCd tt 
h d a d / 1 

i= 1 





(2) 
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Here, Cd = r(d/ 2 +i) ’ 9 ^ ie s P^ n degeneracy factor, V c j t = 2 d nf=i a * is the volume of an d-dimensional rect¬ 
angular whose i-th side has length 2aj. T(Z) = / 0 °° dxx l ~ 1 e~ x is the gamma function and x = j + X)f=i 

The grand potential for the Fermi system, 

q = — ^ Zn(l + zexp(—/3e )) (5) 

€ 

j3 = Tjip, where k being the Boltzmann Constant and z = exp(/3/ii) is the fugacity, where /r being the chemical 
potential. Using the Thomas-Fermi semiclassical approximation[25] and re-writing the previous equation, 


roo 

q = Qo - / p(e)Zn(l + zexp(-/3e )) 
Jo 


So, using the density of states of Eq. (3) we finally get the grand potential, 

q = q 0 + Br(j + l)(kT)Xf x+1 (z) 

where, qo = —ln{ 1 — z) and fi(z) is the Fermi function which is defined as, 


roo p -1 °° 

f P {z ) = / dx 

Jo 2 ^ + 1 ^ J 


jp 


( 6 ) 
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3. Thermodynamics of Fermi gas under generic power law potential in d dimension 

3.1. Number of Particle 
The number of particles N can be obtained, 


da 

N — = N 0 + 


sC„r(f + i)V; nf_i r(i + 1) 


h i h m n <i_ i c j/». 

Here, No = jf— is the ground state occupation number. 

Now, defining, 

d ut i 

Vi = Vdf[( —) 1 /ni r(—+ i) 

±A ; Ci rii 

i=i 1 1 


(KT)Xf x (z) 


A' = 


hbi d/2 + 1 


tt 5 (fcT)i d/l + 1 


]1 /d 


(9) 


( 10 ) 

( 11 ) 


It is noteworthy, 


lim V^Vd 

rii—too 


lim X = T 

n.i—1 oo' l 


lim A / = Aq = 




(2-KmkT) 1 / 2 


( 12 ) 

(13) 
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So, if we choose l = 2 and b = ^ from Eq. (14) we get Ao = (■- >7VlllkT p /2 ? which is the thermal wavelength 

of nonrelativistic free massive fermions. However, it should be noted that, when l ^ 2, X' then depends on 
dimension. With d = 3 and d = 2, thermal wavelength of photons are respectively 0n l/ 2 kT and 
which can be obtained from from Eq. (11) by choosing b = c, where c being the velocity of light. So, 
one can reproduce the thermal wavelength of both massive and massless fermions from the definition of X' 
with more general energy spectrum. But one needs to consider the effects of antiparticles to calculate the 
thermodynamic quantities of ultrarelativistic quantum gasf2?j. 


The number of particle equation is then written as, 

N-N 0 =g^f x (z) (15) 

The number of particle equation for free massive fermions (with l = 2, a = -^- L , all m —> oo) in d 
dimesnsional space can be obtained from Eq. (15), 

N-N 0 = g^f,(z) (16) 

A 0 2 

which gives the exact equation for number of particles at d = 3HM2]. 


3.2. Internal Energy 

From the Grand Canonical Ensemble internal energy E is, 

E = -(§).,K 

- gaj(f + 1) ^ 


= NkT X 


h d b d /l nf=i c i 
fx+ 1(-) 

fx( z ) 


In case of free massive fermions, 


E = NkT 


dfd/ 2+1 (A) 


2 fd/ 2 (z) 

which is in accordance with the exact expression of E for d = smm- 


(17) 

(18) 


(19) 


Now as T —> oo, from Eq. (19) it is seen, the internal energy becomes, E = NkT\■ For free massive 
fermions it is E = ^NkT , which becomes ~NkT , when d = 3. Thus E approaches the classical value at 
high temperature. The exact same trend is also seen in case of Bose gas[l]. 


3.3. Entropy 

The entropy S can be obtained from Grand Canonical Ensemble, 

5 = kT(^) z y-Nklnz + kq 

= Nk ijk(x + l)f x +i(z) - In z\ 


( 20 ) 
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As before, for free massive fermions, one can find Eq. (20) approaches to, 

S = + ~ lnz ] ( 21 ) 
Again at d = 3, Eq. (21) reduces to same expression for entropy as Ref. [TT1 fl2| 


3-4■ Helmholtz Free Energy 

From the Grand Canonical Ensemble we get the expression of Helmholtz Free Energy for Fermi system, 


A = —kTq + NkT\nz 

= -NkT^+lill + NkT\nz 
fx( z ) 

In case of free massive fermions, the above expression reduces like below, 

A U+ii*) , 

NkT ~ fd(z) + nZ 

Now, for d = 3, the above equation produces the exact expression for Helmholtz free Energy [TT1112] . 


( 22 ) 


(23) 


3.5. Pressure 

Rewriting equation (15) stating the number of particles, 


N - N 0 


^nti(f)^r(i- + i) 


N - No g 

= yd9x( z ) 


V'i 


Now a very well known expression for the nonrelativistic d— dimensional ideal free Free gas{TTj, 


N -N 0 
V d 


jdfd/ 2 (~) 


Comparing the above equations, we can say is a more generalized extension of Vd- Where, 

d bT i 

v^ = v d ll (_) 1 Mr(- + i ) 

f 4- Cj n i 

i=i 

It represents the effect of external potential on the performence of trapped fermions. Calling the effective 
volume the grand potential can be rewritten as, 


gVi 

Q = Qo +9-yffx+i( z ) 


(24) 


So, the effective pressure 


p' = 4 (- 8q 


^dvy 


gkT 

\' d 


fx+ 1(«) 


(25) 
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Which can be rewritten as, 


NkTf x+1 (z ) 

K /*(*) 


(26) 


The above equation is very general equation of state for any dimensionality d, any dispersion relation of the 
form (oc p s ) having any form of generic power law trap and obviously it is expected that it will reproduce 
the special case of free system. For free system the equation (26) becomes, 

(27) 

P oV d V d fd(z) 

2 

which is in accordance with Ref. mm at d=3. 


Now, comparing Eq. (18) and (26) we get, 

P'K = f (28) 

A 

For d-dimensional free Fermi gas one can obtain from previous equation 

PV d = -E (29) 

a 

This is an important and familiar relation, PV = |E when d = 3 nn uni m m - This actually shows 
equation (28) is a very significant relation for the Fermi gas irrespective whether they are trapped or free. 
And in case of trapped fermions effective volume and effective pressure plays the same role as volume and 
pressure in current textbooks and literatures. Interestingly, the Bose gas also maintains this equation. |T] 


3.6. Heat Capacity 


Heat capacity at constant volume C v , 


Cy = T(^) N y 

= Nk[x(x + l)^p/ x +1(*) ~ 


(30) 


For Free massive fermions, the expression becomes, 


A A „ A 2 fd(z) 

Cv = i + 


(31) 


And in the high temperature limit of C v approches its classical value as it becomes X-/V& for trapped system 
and ^Nk for free system, which is | Nk, when d = 3. 

Now, heat capacity at constant pressure C p , 

Op = r(||)„,p 

= Nk[(x + 1) 2 /x+i(^)/x-i(^)(^7d) 3 - X(X + 1 )/x+i(^)^7d] 


(32) 
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In case of free massive Fermi gas, the above equation reduces to, 


,d 


d.d 


C P = Nk[(— + 1 YfUW^z)^)* - -(- + l)h +1 (z ) w ] 


2 2 


A ,D 


(33) 


It coincides exactly with Ref. DUE] for d = 3. Again in the high temperature limit C p becomes (x + l)IVfc 
for trapped system and + 1)IVA; for free system, which is | Nk , when d = 3. So, in the high temperature 
limit C p approches its classical value. 

Now, the ratio, 7 = (pA) for Fermi gas is given by, 

i/x+ibd 


7 = 


^ +i ) 2 /x+i( A / r w -x(x+D^f 
x(x + 1 )^t/x+iW - 


(34) 


Now, in high temperature limit, the above equation becomes, 

(x + 1) 2 - x(x + 1) 


7 = 


= 1 + 


1 


X(X + !) - X 2 X 

In case of free system, choosing all n* —> 00 , we get from the above equation, 

l 

7 = 1 + d 

With d = 3 and l = 2, 7 equals |, thus obtaining the classical value at high temperature limit. 

3.7. Isothermal Compressibility 
The Isothermal compressibility of Fermi gas can be obtained, 

K T = ~Vd(j^) N ,T 

= -v'(-) (—) 

v d\ Q z >N,T\Qy/JN.T 

Vd fx- i(z) 


(35) 


(36) 


NkT f x (z) 


(37) 


which reproduces the same result for isothermal compressibility of free massive Fermi gas at d = 3 mi And 
as T —> 00 up takes the classical value for free system, which is -p. 


4. The thermodynamic properties of a degenerate Fermi gas under generic power law 

potential 

At low temperature, can approximate the Fermi function and write it as quickly convergent Sommerfield 
series m 


fp(z) 


(In z) v 

r(p +1) 


[i+p{p- 1) 


7 r 2 1 

6 (lnz) 2 


+ P(P~ l)b-2)(p-3) 


7vr 4 1 
360 (In z) 4 


+ ...] 


(38) 
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At T = OK, taking only the first the first term in Eq (38). Substituting this into Eq. (15) we get 

g c n T(i+ i)v d nUn± + i) 


N — Nq = N r = 


Which turns out, 


Ep = 


h d b d / 1 ntic! /n T( X + l) 

h d b d / 1 nti4 Kr (x+i)iv e i 


-E X F 


(39) 


gC n r(% + i)v d nf = i r(i + i) 

Following the method of Ref. mmm we approximate the chemical potential and fugacity from Eq.(15) 

7r 2 . kT 


(40) 


h = fcTlnz = E F [1 - (x - l)^)^) 2 ] 


(41) 


Using these approximation we can calculate the thermodynamic quantities of the previous section, 

X 


E 

N 

S_ 

Nk 


X + 1 6 Ep 

9 

X7T 


— = 777 ~kT 

3 Ep 
EpN 


P = 

Cv_ 

Nk " 

Kp = 


r TT 2 .kT o 

(x + i)v [ +(x + 1) T ( e7 )1 


■fcT 


_ X7T 
3Ep 


In case of free massive fermions (choosing 1 = 2), Eq. (42)-(46) becomes, 

d 


E _ 
N ~ 
S 

~Nk 
P = 


„ r ,d , 7r 2 , kT . 9l 
—^ F [!+ (- + !)-(—) ] 


■kT 


dp 2 
6 Ep 

2EpN d ,.p 2 .kT o 

(Wr [1 + ( 2 + 1 »I^ 

A = 

IV/c 6 Ef 

U'd d p 2 kT 2 

KT =mE-p {1 + ^-2^E-p )} 


(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


At temperature T = OK, entropy S = 0 which is accordance will 3rd law of Thermodynamics. The internal 
energy, pressure and isothermal compressibility T = OK, 


E 0 

Po 

kto 


X 

X + 1 




1 N 

(x+Tiu 7 


Ep 


^x 


iye f 


(52) 

(53) 


(54) 
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In case of Free massive Fermi gas the above equations reduces to, 


Eq = 

d NE f 
d + 2 

(55) 


2 N 


Po = 

Ef 

(d + 2)V F 

(56) 


V d 


Kto = 

NE f 2 

(57) 


At d = 3, Eq. (55) and (56) become exactly same as in Ref. [IT] 

5. Discussion 

Thermodynamics of ideal Fermi gas in the presence of an external generic power law potential are dis¬ 
cussed in this section. It is seen the effective volume V' d is a very salient feature of trapped system, playing 
the same role in trapped system as the volume in free system, which enables us to treat trapped fermi gas 
as well the bose gasp.j to be treated as a free one. Difference between Fj and Vd is that, the former depends 
on temperature and power law exponent while the latter does not. But as all ni —> oo, approaches 
Vd . In this process the more general thermal wavelength X' is defined with arbitrary kinematic parameter 
in any dimension. It was shown how X' can reproduce the thermal wavelengths of literatures in different 
dimensions. In case of trapped Fermi gases, and A' enable all the thermodynamic functions of the system 
to be expressed in a compact form similar to those of free Fermi gas. 

At first the density of states and grand potential is calculated in section 2 . All the thermodynamic quantities 
are derived from the grand potential in section 3. It is seen that all the thermodynamic quantities for trapped 
Fermi gas are some function of Fermi functions, just as in the case of free Fermi gas. But in the former case 
the Fermi functions depend on X = j + TT:=i and z, whether in the later case Fermi functions depend on 

j and z. And as n* —> oo, the mathematical form of thermodynamic quantities of trapped system reduce 
to that of free system. Nevertheless it is noteworthy that Eq. (28) is a very remarkable relation for quantum 
gases as both Bose[T] and Fermi system maintains it. 


In general, the thermodynamic quantities of trapped system differ from free system. We can specifically 

check this by comparing free system with harmonically trapped potential. Let, d = 3, a = Z = 2, 

2 

n\ = u -2 = U 3 = 2, a = and g = 2. Results of some of the physical quantities have been listed in Table 
1. From the table it is seen, thermodynamic quantities are affected by the trapping potential. The signa¬ 
ture of trapping potential is present in the low as well as in the high temperature limit of thermodynamic 
functions. As we know chemical potential approaches Fermi Energy as T —> 0 in case of free Fermi gas, 
same phenomena is also observed in case of trapped Fermi gas, although Fermi Energy do vary comparing 
trapped system with the free one. Let us turn our attention to low temperature limit of Fermi gas. It is 
seen both Cy and S has same numerical value in this limit just like the free system and goes to zero at 
T = OK . The later actually is a manifestation third law of thermodynamics. But most significantly internal 
energy and pressure of trapped Fermi gas do not go to zero as T = OK just as free Fermi gas. According to 
equation (46), no matter what power law exponent one chooses, ground state pressure never goes to zero. It 
suggests, the ground state energy and ground state pressure seen here is clearly a quantum effect arising due 
to Pauli exclusion principle due to which the system can not settle down into a single energy state as in the 
case of Bose gas. And therefore spread over a lowest available energy states. More interestingly, isothermal 
compressibility of Fermi system is nonzero at T = OK, indicating zero point pressure is not merely a constant 
but depends on volume. 
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Table 1. Dissimilarity between free and harmonic-potential-trapped nonrelativistic Fermi gases in d = 3. 


Physical quantity 

Free gas 

Trapped gas 

Fermi Energy 

fAAfh) 

hu:{2,N) 1 ^ 

Fermi Temperature 


^(3 A) 1 / 3 

Internal energy 

, / 5 (z) 

l NkT fl( Z ) 

3iVMtg 

Internal energy at lower temperature 

fACMl + f^g) 2 ) 

|JVE F (1 + 3|1(|C)2) 

Ground state energy 

f NEp 

f NE f 

Ground state pressure 

2 N rp 

5 V^ F 

1 N p 

4 V'^ F 

Particle number at ground state 

J^(2 mE F f/ 2 
|NkT 

1 / E f \3 

3 1 hoj ) 

Internal energy at higher temperature 

3 NkT 

Isothermal Compressibility 

V / 1 / 2 O) 

NkT 73 / 2 ( 2 ) 

V' h(z) 

NkT f 3 (z) 

Isothermal Compressibility at higher 
temperature 

1 

P 

1 

P' 

Isothermal Compressibility at lower tem¬ 
perature 

3V n 7T 2 ! kT \ 2l 

2 NE F P 12\E F ) i 

3V n tt 2 ( kT ^1 

NEp P 3 V E F > \ 

Isothermal Compressibility at T = OK 

3 V 

2 NEp 

3V 

NEp 

Chemical potential at lower temperature 

Ef(1- fj(f) 2 ) 

m i-T(f) 2 ) 

Chemical potential at T = OK 

Ep 

Ep 

Chemical potential at higher temperature 


-fcTln[6(g) 3 ]) 

C'v 

(8 /C> 

' 4 ) f 3 (z) u//i (z) 

'^ ’’J 

1<l/4(z) n/3(z) 

1 Z / 3 (z) J / 2 (z) 

Nk 

jtj: at lower temperature 

7T 2 kT 

2 Ep 

vr 2 AT 

Ep 

at higher temperature 

3 

_2_ 

3 


6. Conclusion 

From the grand potential, the thermodynamic properties of Fermi gas trapped under generic power 
law potential have been evaluated. The calculated physical quantities reduce to expressions available in 
the literature, with appropiate choice of power law exponents and dimensionality. The thermodynamic 
quantities are studied further closely in the degeneracy limit. It is found pressure, energy and isothermal 
compressibility is nonzero, with any trapping potential, indicating the governing power of Pauli exclusion 
principle. In this manuscript, we have restricted our discussion in case of ideal system under trapping 
potential. It will be vey interesting to see the effect of interaction in the degeneracy limit. 
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